Algebraic fidelity decay for local perturbations 
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From a reflection measurement in a rectangular microwave billiard with randomly distributed 
scatterers the scattering and the ordinary fidelity was studied. The position of one of the scatterers 
is the perturbation parameter. Such perturbations can be considered as local since wave functions 
are influenced only locally, in contrast to, e.g., the situation where the fidelity decay is caused by 
the shift of one billiard wall. Using the random-plane-wave conjecture, an analytic expression for 
the fidelity decay due to the shift of one scatterer has been obtained, yielding an algebraic 1/t decay 
for long times. A perfect agreement between experiment and theory has been found, including a 
predicted scaling behavior concerning the dependence of the fidelity decay on the shift distance. The 
only free parameter has been determined independently from the variance of the level velocities. 

PACS numbers: 05.45.Mt, 03.65.Nk, 03.65.Yz 



A versatile tool to describe the quantum mechanical 
stability of a system against perturbations is fidelity in- 
troduced by Peres 25 years ago The fidelity ampli- 
tude is defined as the overlap integral of an initial state 
with itself, after having experienced two slightly dif- 
ferent time evolutions, 



f{t) = {Me 



IV-o). 



(1) 



Here Hq and Hi — Hq + XV are the two relevant Hamil- 
tonians, and parameter A describes the strength of the 
perturbation. Throughout this letter the time will be 
given in units of the Heisenberg time tn = h/AE, where 
AE is the mean level spacing. The fidelity F{t) is the 
modulus square of the fidelity amplitude. The concept 
of fidelity is closely related to the concept of decoher- 
ence which is widely used in the community of quantum 
computation [1, 

For very weak perturbation strengths A, in the pertur- 
bative regime, Eq. H]) reduces to 



fit) = (V'ol 



IV-o) 



(2) 



where Vdiag is the diagonal part of the perturbation in 
the basis of Hq. Assuming a Gaussian distribution for 
the matrix elements of Vdiag Eq- @ describes a Gaussian 
decay of the fidelity beyond the Heisenberg time. With 
increasing perturbation strength the Gaussian decay is 
taken over by an exponential decay below the Heisen- 
berg time with a decay constant which can be calculated 
from the matrix elements of the perturbation by means 
of Fermi's golden rule [1, @, @. For perturbations with 
missing diagonals the Gaussian decay is suppressed leav- 
ing only the exponential decay, a situation termed freeze 
0. For very strong perturbation strengths the fidelity 
decay becomes independent on the perturbation strength 
and reflects the classical dynamics. In chaotic systems an 
exponential decay with the classical Lyapunov exponent 
is observed Q, whereas in integrable and diffusive sys- 
tems an algebraic tail is found [3, • For more details 
it is referred to the review article by Gorin et al 



In all references mentioned above global perturbations 
have been considered, meaning that there is a total re- 
arrangement of spectrum and eigenfunctions already for 
moderate perturbation strengths. For this situation ran- 
dom matrix theory, assuming uncorrelated Gaussian dis- 
tributions for the matrix elements of V, has yielded a 
number of perturbative and exact results 13, hj. The fi- 
delity decay by global perturbations has been studied ex- 
perimentally in a microwave billiard, where one wall was 
shifted [3], and in a vibrating aluminum block, where 
the temperature took the role of the perturbation pa- 
rameter [H, 16|. In chaotic systems the long time fidelity 
for global perturbation always is either Gaussian or sin- 
gle exponential, a very unsatisfactory situation, e.g., for 
quantum computing. 

The role of local perturbations, realized in billiard sys- 
tems, e.g., by the shift of an impurity, has been more 
or less ignored in the past. This is somewhat surprising, 
since many, if not most of the perturbing interactions in 
real systems are short ranged. Examples are diffusive 
jumps or flips of neighbored spins in solids leading to 
the decay of s pin echoes in nuclear magnetic resonance 



(see e.g. Ref. [17[). Another example is the twinkling of 



stars, caused by the diffusive motion of the atoms in the 
atmosphere. 

In a previous work we studied the spectral level dy- 
namics in microwave billiards for global variations, real- 
ized by a shift of one wall, and local variations, realized 
by the shift of an impurity [l^. For global variations a 
Gaussian eigenvalue velocity distribution was found, as 
expected for chaotic systems. For local perturbations the 
situation was found to be completely different. The in- 
troduction of a point- like impurity leads to a shift of the 
eigenvalues by 



AEn ^ a\ilJn{r)\ 



(3) 



Hi. 



where a describes the strength of the perturbation, and 
^pn is the wave function at the scatterer position in the ab- 
sence of the scatterer. Using the random-plane- wave as- 
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FIG. 1: Part of the reflection spectrum for two perturber 
positions differing by |Ar| = 4 mm. The inset shows a sketch 
of the used microwave biUiard (scatterers are shown three 
times larger). The antenna position is marked by a cross. 



sumption a modified Bessel function was obtained for the 
velocity distribution. With increasing radius of the im- 
purity or, equivalently, with increasing frequency a grad- 
ual transition from local to global behavior was found. 
Somewhat later our results had been verified in super- 
symmetry calculations [l9| . 

In view of the fact that the velocity distributions for 
global and local perturbations are completely different 
one could expect that the same is true for the fidelity 
decay. It will be shown in this letter that this is really 
the case. We shall see that for local perturbations the 
fidelity decays algebraically for long times, i. e. the situ- 
ation is much more favorable than it is the case for global 
perturbations. 

The fidelity in its original definition ([T]) is hardly ac- 
cessible experimentally. This was our motivation in a 
previous work fl3| to introduce the scattering fidelity be- 
ing defined as 



(4) 



Here Sab is a scattering matrix element, and 
C[SabT S\^*]{t) the Fourier transform of the cross- 
correlation function of the respective matrix elements. 
It was stated in Ref. that for chaotic systems and 



weak coupling of the measuring antenna the scattering 
fidelity approaches the ordinary fidelity. In the present 
work it has been possible to study both quantities at the 
same time and to compare the results. Since in all previ- 
ous experiments on fidelities actually scattering fidelities 
or related quantities had been measured, though without 
stating this explicitly, such a comparison is of fundamen- 
tal interest. 

The system studied was a rectangular microwave bil- 
liard with height h = 8 mm, side lengths a = 340 mm. 



b = 240 mm and 19 disks with a diameter oi D = 4.6 mm 
placed randomly inside the billiard (see inset of Fig. [T]). 
Another disk of the same size has been moved in steps 
of |Ar| = 1mm through the billiard. With a fixed wire 
antenna we measured the reflection spectrum for 300 dif- 
ferent positions of the moving disk in a frequency range 
from 3.5 GHz to 6 GHz. In this frequency range the bil- 
liard can still be treated as two-dimensionally, and there 
is a complete agreement with the corresponding quan- 
tum mechanical system [ici] . Figure [T] shows part of the 
reflection spectrum for two slightly different positions of 
the movable disk. In addition we measured the eigenfunc- 
tions [IJ] to make sure that they are delocalized and their 
intensities Porter-Thomas distributed. Thus the system 
is fully chaotic and ballistic in the investigated frequency 
range. In the measured frequency regime it was also pos- 
sible to extract resonance positions and amplitudes from 
the spectrum by a Lorentz flt. These quantities will be 
needed later for the determination of the ordinary fidelity. 

First we shall derive an explicit expression of the fi- 
delity decay caused by the shift of a local perturber. We 
assume that the system is completely chaotic (in the ex- 
periment this was verified by measuring the eigenfunc- 
tions, see above). In this case we may average expression 
([1]) over all possible initial states resulting in 



/W = ^Tr(e 



2TTlHot~2TTlHit\ 



(5) 



where N is the number of states taken in the trace. In 
the present case Hq and Hi correspond to the Hamil- 
tonian of the billiard with the perturber placed at two 
different positions. For a weak, point-like perturbation, 
the perturber just produces a small shift of the eigenen- 
ergy proportional to the intensity of the unperturbed 
wave function at the perturber position, see Eq. ([3]), i. e. 
the Hamiltonian in the basis of the unperturbed system 
is given by 

Hnmir) = 6nrniE^ + a\Mr)\'') , (6) 

where the are the eigenenergies of the unperturbed 
system. For this approach to be valid it is essential that 
the shift induced by the perturber is always small com- 
pared to the mean level spacing, i. e. we never leave the 
perturbative regime. We then obtain from Eqs. ^ and 
® 



fit) = ( 



_ /g277ia(|Vir 



(7) 



where the abbreviation tpi = tp{ri) has been used. The 
brackets denotes the average over all N participating 
states. To calculate this average, we now apply the ran- 
dom plane wave assumption. The average is most conve- 
niently expressed in the following way [22| , 



fit) 
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(8) 
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FIG. 2: Scattering fidelity (color online) for difi^erent shifts 
of the perturber |Ar| = 1mm (open circles, black), 2mm 
(crosses, blue), and 4mm (filled squares, red). The solid line 
corresponds to the theoretical prediction of Eq. (|10p . The 
results were obtained for the frequency range from 3.5 GHz 
to 6 GHz. The insert shows the same results in a double- 
logarithmic plot. 

where K is a (2 x 2)-niatrix of which the inverse 

K-^ = ( ^^^^A (9) 

can be expressed in terms of two-point correlation func- 
tions only. Using again the model of random plane waves 
[23II the two-point correlation function can be calculated 
yielding {ipiipj) — ^ Jo(fc|ri — |), where A is the billiard 
area. Performing the integrations one obtains the final 
expression for the fidelity amplitude 

/(t)= [l + (At)^]"^ , (10) 

where 

X = ^^l^j2(^k\Ar\). (11) 

|Ar| is the shift of the scatterer. For large t the fi- 
delity amplitude decays algebraically with f{t) ^ 1/t. 
The only free parameter in Eq. (jlip is a. This parame- 
ter can be obtained independently from the variance of 
the level velocities (see Ref. d). This allows us to com- 
pare the experimental results with theory without any 
free parameter. 

Figure [5] shows the scattering fidelity, as obtained 
from Eq. ([4]), for three different perturber shifts. The 
solid lines correspond to the theoretical prediction from 
Eq. (llOp . A good agreement between theory and exper- 
iment is found in all cases, apart from minor systematic 
deviations for small times. Values obtained for the cou- 
pling parameter a by a fit deviate from those from the 




FIG. 3: Same as Fig. (2] but plotted on a rescaled time axis 
using the scaling property of Eq. (fTOjl . 

level velocities by a few %. The inset shows the results in 
a double-logarithmic plot illustrating the algebraic decay 
for long times t. Equation PH)) exhibits a scaling behav- 
ior: On a rescaled time axis Xt all experimental results 
should fall onto one single curve. Figure [3] demonstrates 
that this is indeed the case. 

So far we have discussed the results for the scattering 
fidelity. Let us now see, how the ordinary fidelity H]) can 
be obtained. Since the measurement has been performed 
at a fixed antenna position tq, the initial state is localized 
at ro, i.e. \->po) = ko) = ^ir - ro). Expanding IV'o) in 
terms of eigenfunctions {ipn) of -^ii ^^d {^po\ in terms of 
eigenfunctions (t/;"! of Hq, Eq. ^ reads 

fit) = 5:^°(ro)^i,(ro)e^-<^-<)*(^5^|^^) . (12) 

n,m 

The shift of the energy is a first order effect of the per- 
turber, the change of the eigenfunctions being of the next 
order. Neglecting these higher order effects, we may ap- 
proximate (V'nIV'm) ~ ^nm tO obtain 

/W-E^n('^0)'^n(-0)e^"(^"-^")'. (13) 
n 

All quantities entering the sum on the right hand side 
are available from the experiment, the eigenenergies from 
the resonance positions, and the eigenfunctions at the 
antenna position from the resonance depths |i24]. Al- 
together 64 resonances have been taken. In Fig. |4] the 
results for the fidelity are presented. A very good agree- 
ment between theory and experiment is found. Here the 
relative deviation for the coupling parameter a taken 
from the variance of the level velocities and the one ob- 
tained from the fidelity even agree up to 1%. No smooth- 
ing has been applied. The ordinary fidelity, too, shows 
the scaling behavior predicted by Eq. (jTU]) as is pro- 
nounced in Fig. [51 In addition the collected results for 
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FIG. 4: Ordinary fidelity determined as obtained from the 
resonance position and depths, see text. The meaning of the 
symbols is the same as in Fig. (2] 
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FIG. 5: Scattering fidelity, obtained from a superposition of 
the results from all perturber shifts (blue crosses) and ordi- 
nary fidelity (red squares) on a rescaled time axis. The solid 
line corresponds to the theoretical prediction. 

the scattering fidelity, obtained by an averaging over all 
perturber shifts, is plotted, showing the quality of the 
agreement between both types of fidelity. 

In summary we have shown that for a local perturba- 
tion, realized by the shift of a perturber in a microwave 
billiard, the fidelity decays algebraically with a.l/t long- 
time behavior. This is in contrast to the exponential 
or Gaussian long-time behavior observed in the fidelity 
decay due to global perturbations. All results could be 
quantitatively explained within the random plane wave 
model, including a scaling prediction on the dependence 
of the fidelity decay on the perturber shift. In addition we 



could show that scattering fidelity and ordinary fidelity 
are identical within the experimental errors. 
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